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ABSTRACT 

Gradual shear acceleration of energetic particles in gamma-ray burst (GRB) 
jets is considered. Special emphasis is given to the analysis of universal structured 
jets, and characteristic acceleration timescales are determined for a power-law and 
a Gaussian evolution of the bulk flow Lorentz factor 7^ with angle <fi from the 
jet axis. The results suggest that local power-law particle distributions may be 
generated and that higher energy particles are generally concentrated closer to the 
jet axis. Taking several constraints into account we show that efficient electron 
acceleration in gradual shear flows, with maximum particle energy successively 
decreasing with time, may be possible on scales larger than r ~ 10 15 cm, provided 
the jet magnetic field becomes sufficiently weak and/or decreases rapidly enough 
with distance, while efficient acceleration of protons to ultra-high energies > 10 20 
eV may be possible under a wide range of conditions. 

Subject headings: gamma rays: bursts — acceleration of particles 



1. Introduction 

There is mounting evidence today that gamma-ray bursts (GRBs) are associated with 
collimated relativistic outflows or jets (Rhoads 1999; Kulkarni et al. 1999; Greiner et al. 
2003). It is still a matter of ongoing research, however, what types of internal structures are 
actualized within these jets. Within the fireball framework, for example, two different kinds 
of jet models seem to be compatible with the observations, i.e., (i) the uniform "top-hat" jet 
model (Rhoads 1999; Frail et al. 2001; Lamb et al 2004), in which all the hydrodynamical 
quantities (e.g., Lorentz factor, energy density) are essentially the same within some well- 
defined opening angle 4>j around the jet axis but sharply drop outside of <fij, and (ii) the 
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universal structured ("power-law" or "Gaussian") jet model (Rossi et al. 2002; Zhang & 
Meszaros 2002; Kumar & Granot 2003; Zhang et al. 2004), in which the hydrodynamical 
quantities are rather smooth functions of the angle <fi from the jet axis. 
Here we analyze the implications of such jet structures for the acceleration of energetic 
particles. We focus on shear acceleration as a promising mechanism for converting the kinetic 
energy of the flow into nonthermal particles and radiation. Such a mechanism has previously 
been successfully applied to the relativistic jets in Active Galactic Nuclei (Ostrowski 1998, 
2000; Rieger & Duffy 2004, 2005a,b). Shear acceleration is based on the simple fact that 
particles may gain energy by scattering off (small-scale) magnetic field irregularities with 
different local velocities due to being systematically embedded in a collisionless shear flow (cf. 
Rieger & Duffy 2005b, for a recent review). In the case of nonrelativistic shear acceleration 
it is straightforward to show that local power-law particle momentum distributions f(p) oc 
p-(3+a) can k e g enera ted, assuming a momentum-dependent mean scattering time of the 
form r oc p a with a > (cf. Berezhko & Krymskii 1981; Rieger & Duffy 2005a). Hence, 
for r scaling with the gyro-radius, i.e., r oc p, this results in a power-law particle number 
density n(p) oc p 2 f(p) oc p~ 2 , and thus a synchrotron emissivity j v oc z/ -1 / 2 . 



2. Shear acceleration in structured GRB-type jets 

Let us consider an idealized radial, relativistic flow profile, which in four vector notation 
is given by (cf. Rieger & Duffy 2005c, for the more general case) 

u a = lb (M r (0)/c,O,O) , (1) 

where a — 0, 1, 2, 3, and <fi denotes the polar angle in spherical coordinates, and 7b = 7&(0) = 
[1 — f r (0) 2 /c 2 ]~ 1 / 2 is the bulk Lorentz factor of the flow. In the comoving frame, the related 
gradual shear acceleration coefficient can be cast into the form (e.g., Webb 1989, eq. 3.27) 

where p' denotes the comoving particle momentum, t' ~ A'/c is the mean scattering time 
and r is the relativistic shear coefficient. We are interested in the strong scattering limit 
(corresponding to oo g ' r' <C 1, with uo' g being the relativistic gyro-frequency measured in the 
comoving jet frame, i.e., assuming a sufficiently weak longitudinal mean magnetic field and 
the presence of strong turbulence, so that collisions are efficient enough to restore isotropy), 
where the shear coefficient is given by V = (c 2 /30) <j a p a af3 (see Webb 1989, eq. 3.34) and a a p, 
with a, p — 0, 1, 2, 3, is the usual covariant fluid shear tensor (cf. Webb 1989; Rieger & Duffy 
2004, 2005c, for more details). Using spherical coordinates and the velocity profile above, it 
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can be shown that the relativistic shear coefficient becomes (cf. Rieger & Duffy 2005c) 



vl 3 2 ( dv r 

^2 + 4^2 ^ 



(3) 



which in the nonrelativistic limit ( 7 & — > 1) reduces to the (nonrelativistic) viscous transfer 
coefficient derived by Earl et al. (1988) (their eq. 7) when the latter is expressed in spherical 
coordinates and the corresponding velocity profile v = v r (<p) e r is applied. The (comoving) 
timescale t acc ~< p' > / < p' > for the shear flow acceleration of particles then becomes 

45 c r 2 
tacc (r ' 0) " 4(4 + a) A 7 7 2 k 2 + O.75 7b 2 (<%r/<90) 2 ] ' ( } 

where r is the radial coordinate measured in the cosmological rest frame (i.e., the supernova, 
collapsar or merger rest frame), and where a power- law dependence r' = A'/c oc p' a has been 
assumed. As the jet flow is diverging, t acc obviously increases with r to the square. 
In order to investigate the acceleration potential of structured GRB jets, two applications 
appear particularly interesting (e.g., Zhang & Meszaros 2002; Kumar & Granot 2003; Zhang 
et al. 2004): a power-law model, in which 7^ is a power-law function of <p outside a core of 
opening angle C , i.e., ^ h {<j>) = 1 + (7&0 - 1)(1 + [0/0 c ] 2 )" 6/2 , with 1.5 < 6 ~ 2 (cf. Zhang 
& Meszaros 2002), and a Gaussian model with 7 &(0) = 1 + (7&0 — 1) exp[— 2 /20 2 ], where 
7m denotes the Lorentz factor at the jet axis, and typically <p c = 0.1 rad (Zhang et al. 
2004). The shear acceleration timescale (with the time and space coordinates measured in 
the cosmological rest frame and particle Lorentz factor measured in the comoving jet frame) 
then becomes 

, , r M = -J^- j (§ + ! Jf&^mJw g) " Power - law 

" 4(4+a) c ***r 1 (sM^fcMr Gaussian model ' 

(5) 

Note that in contrast to nonrelativistic parallel shock acceleration, t acc oc 1/A as the probabil- 
ity of a particle sampling a higher shear, and thus a more energetic scattering event increases 
with A. The general evolution of t acc /t as a function of 0, with t = 45r 2 /[4(4 + a) A'c], is 
illustrated in Fig. 1, indicating that a power-law model becomes somewhat more favourable 
than a Gaussian model. Obviously, the shear acceleration timescale in structured jets gen- 
erally increases with <fi due to the decrease in 7^, suggesting different maximum attainable 
energies in the local comoving frame. Hence, at constant r the higher energy particles and 
thus the higher energy emission are expected to be naturally concentrated closer to the jet 
axis, i.e., toward smaller (provided the magnetic field is not very inhomogeneous across the 
jet), enforcing the effect that observers looking at the same GRB from different directions 
will see different gamma-ray light curves. 
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In general several constraints need to be satisfied in order to allow for efficient particle 
acceleration in GRB jets, e.g., the particle acceleration process may be limited by (i) radiative 
synchrotron losses, (ii) the transversal confinement size, (iii) particle escape via cross-field 
diffusion or (iv) the duration of the jet expansion. The first constraint implies that the 
acceleration timescale has to be smaller than the comoving synchrotron cooling timescale, 
which, for particles with Lorentz factors 7' and isotropic pitch angle distribution, is given by 

9m 3 c 5 f , , 9 /1\ /m\ 3 /10 3 G\ 2 M 

with m p being the proton mass. According to eq. (5) the minimum shear acceleration 
timescale is 

5 r 2 

tacc ~ 2 ^1 ' (7) 

so that for a gyro-dependent particle mean free path A' = £ r' , with £ < f, (as motivated 
by the observationally implied particle spectra with momentum index ~ —2) the ratio of 
shear to cooling timescale becomes independent of the particle Lorentz factor. Note that 
eq. (7) is also the timescale for a uniform flow with v r independent of (cf. eq. [4])! For a 
given magnetic field strength the acceleration process may thus work efficiently as long as 
the velocity shear remains sufficiently high, i.e., as long as the radial coordinate satisfies 

suggesting that, in contrast to the acceleration of protons, efficient electron acceleration is 
suppressed in the presence of high magnetic fields. It seems very likely, however, that due to 
the expansion of the wind, the comoving magnetic field B' will depend inversely on r, i.e., 
B' oc 1/r" with f3 > 0. In order to study possible implications in more details, we may thus 
consider the following simple parametrization for the commoving jet magnetic field strength, 
B' = 1 000 b (1 13 cm/r) 13 G. For (3 = 1 and b = 30 this expression corresponds to the lower 
limit required by Waxman (1995) in order to allow for efficient (second order) Fermi-type 
acceleration of protons to ultra-high energies during expansion of the wind (cf. his eq. [4a]). 
While it can be shown then that proton acceleration is nearly unconstrained by radiative 
synchrotron losses, electron acceleration can only work efficiently if the field decays rapidly 
enough and becomes comparatively weak, e.g., for (3 = 2 efficient electron acceleration may 
occur for r > lO 15 ^ 1 / 4 &o /4 (300/7m)) 1/2 cm - 

The second constraint requires the particle mean free path A' to be smaller than the transver- 
sal width it! ~ r(f)j of the jet (with <f>j > <p c being the jet opening angle in the cosmological 
rest frame). For A' = £ r' , £ < 1 and r' g = ^'mc 2 / \e B'), this implies an upper limit for the 
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maximum possible (comoving) particle Lorentz factor given by 

Wx u vm ; ^i 3 g; Vl0 13 cmy V0-3rad; K ) 

Note again, that for B' = 1000 b (10 13 cm/r)^ G, this translates into maximum particle 
energies of E ~ 2.8 • 10 20 £ _1 6 (10 13 cm/rf- 1 (0^/0.3 rad) (7 fe /300) eV in the observer frame. 
For a simple linear scaling, i.e., /3 — 1, the maximum energy would become independent of 
r, while for f3 > 1 the higher energy particles would originate on smaller scales. 
The third constraint requires that the timescale for (cross-field) diffusion in the comoving 
frame t c f ~ R 2 /k', with k' ~ c A'/3, is larger than t acc , i.e., t acc < 3 (r 2 2 )/(c A'). Comparison 
with eq. (5) reveals that t acc scales with r and A' in the same way as t c f, indicating that the 
third constraint may be easily satisfied as long as 7&(0) is larger than a few, cf. Fig. (1). 
Finally, the fourth constraint requires t acc to be smaller than the time t' w = A'/c needed to 
transverse the overall radial (comoving) width A' of the flow, which translates into a lower 
limit for the required particle Lorentz factors. This width is essentially determined by the 
activity of the inner engine driving the GRB and, when specified in the observer frame, should 
be at least of order of the observed GRB duration t A (cf- Piran 2005), i.e., A = t a c = A'/7b- 
Allowing for expansion effects on larger scales with A' ~ r/^y b (cf. Meszaros & Rees 1993), 
we have t' w = max{7btA, r /[lbc]}. Using Eq. (7) one thus finds 

^ - «r e) die) (v) 3 do^) 2 * - { S=) M • 

where x( r ) — 2.7- 10 3 (7fe/300) 4 (10 13 cm/r). For the chosen magnetic field parametriza- 
tion B' = 1000 6 (10 13 cm/r) /3 G, this results in 7 min ~ 10 4 ^(rrip/m) b (r/10 13 cm) 2 "' 3 
(300/7&) 3 x min{(10 sec/t^), x(r)}, indicating that for (3 = 2 the leading term becomes inde- 
pendent of r, which may be favourable for efficient electron acceleration. A more restrictive 
condition, roughly corresponding to 7' nin as implied by x( r ) m ec i- (10)) i s usually associ- 
ated with the requirement to overcome adiabatic losses. Accordingly, shear acceleration can 
only act efficiently on particles that were sufficiently accelerated through other processes. 
In the GRB context it is possible that such particles can be provided, for example, by the 
mechanism responsible for the prompt burst of emission (e.g., shock acceleration). 

In summary, if the jet magnetic field is sufficiently weak (e.g., b ~ 1) and/or decays 
rapidly enough (say, e.g., with (3 — 2), both the first and the fourth constraint may be 
satisfied even for electrons, so that the maximum energy is essentially determined by the 
third constraint. Particles may then, for example, be accelerated efficiently at distances larger 
than r ~ 10 12 (m p /m) cm (cf. constraint [i]), with possible maximum energies, measured 
in the observer frame, of less than ~ 10 18 eV for electrons and ~ 10 21 eV for protons, 
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and successively (i.e., linearly for (3 = 2) decreasing with time. For a high (e.g., b ~ 10) 
and slowly decaying (say, e.g., (5 < 1) jet magnetic field, on the other hand, efficient shear 
acceleration of electrons is virtually excluded, while protons again may well be accelerated 
up to energies > 10 20 eV. 



3. Comparison with shock acceleration 

It is widely believed that shock-accelerated electrons are responsible for the observed 
prompt GRB and afterglow emission via synchrotron radiation processes (e.g., Piran 2005). 
In particular, diffusive electron acceleration at midly relativistic internal shocks (with Lorentz 
factor T s ~ a few) , caused by velocity variations in the relativistic outflow is usually thought 
to be behind the powerful burst of 7-rays (Rees & Meszaros 1994). In the case of nonrel- 
ativistic shocks Fermi acceleration leads to power-law particle spectra N(pj) oc 7 _s , which 
are only dependent on the shock compression ratio p = U\ju 2 (where 1 < p < 4), i.e., 
s — (p + 2)/(p — 1), so that for strong shocks (p = 4) the famous s = 2 result is obtained 
(e.g., Drury 1983; Kirk & Duffy 1999). In general, the acceleration timescale for diffusive 
shock acceleration depends on both the upstream and downstream residence times. For an 
unmodified nonrelativistic shock, one thus obtains (e.g., Drury 1983; Jokipii 1987) 

3 (ki _ _3_p (ki+pk 2 ) , * 

u 1 -u 2 \u 1 u 2 J (p-1) uj 

with ui and u 2 the upstream and downstream flow velocities measured in the shock frame, 
respectively, K\ and k 2 the upstream and downstream (spatial) diffusion coefficients, and u s 
the shock speed as measured in the upstream frame. If the acceleration process operates at 
nearly the Bohm limit [i.e., k ~ Ac/3 = r]r g c/3, with rj ~ 0(1)], one finds 

£acc~6?7r 9 — , (12) 
uj 

assuming K\ ~ k 2 . When u s ~ c this is comparable to the Larmor time which is a result 
also predicted by the simulations of Lemoine & Pelletier (2003). It has been often suggested 
that the nonrelativistic limit can also be approximately applied to midly relativistic shocks in 
GRBs, proposing, for example, that efficient proton acceleration to ultra-high energies might 
be possible (e.g., Waxman 2004). Equating the acceleration timescale (eq. [12]), replacing 
r g by r' g and u s by u' s as the (comoving) internal shock speed, with the cooling timescale 
eq. (6), using B' = 1000 b (10 13 cm/r) /3 G and u' s = (3 s c with (3 S < 1, gives 



' max 
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Note that magnetic field amplification in the vicinity of a shock may yield a field value bo 
well above the usual (background) flow magnetic field, thus leading to a somewhat more 
restrictive condition. In addition, the acceleration timescale also has to be smaller than the 
comoving time ~ r/c^ needed to transverse the width of a shell, which gives 

A more detailed comparison of shear with shock acceleration may perhaps be reached by 
analyzing a critical particle Lorentz factor 7^ oc b^r 1 " 13 (m p /m), defined by equating eq. (7) 
with eq. (12), above which shear acceleration may become more efficient than shock accel- 
eration. Several applications with respect to the acceleration of protons are illustrated in 
Fig. (2), suggesting the possibility that under a reasonable range of conditions shear accel- 
eration in GRB jets may become more relevant for the production of UHE cosmic rays than 
shock-type acceleration processes. Concerning electron acceleration in weak magnetic fields 
(e.g., bo < 1, (3 = 2, r ~ 10 15 cm), on the other hand, shock-type processes would allow for a 
much quicker particle energization (but would also be more severely limited by eq. [13] and 
eq. [14]) and shear effects would only become dominant above 7' ~ 10 8 (10 15 cm/r). 



4. Conclusions 

Using an idealized model we have analyzed the possible role of shear acceleration for 
the energization of particles in relativistic GRB-type jets. Our results suggest that efficient 
electron acceleration on scales r ~ 10 15 cm, with maximum energy decreasing with distance, 
may be possible in the presence of weak magnetic fields, assuming that high energy seed 
particles are provided by the mechanism responsible for the prompt burst of emission. This 
may result in a weak and long duration component in the GRB emission. Protons, on the 
other hand, may reach UHE energies > 10 20 eV under a broad range of conditions. 

We are grateful to R. Blandford for discussions and the referee for very useful comments. 
FMR acknowledges support through a Marie-Curie Individual and a Cosmogrid Fellowship. 
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Fig. 1. — Left: Illustration of the evolution of the bulk Lorentz factor 7& with angle <\> for a 
power-law (b = 1.8) and Gaussian model, respectively, using 7 &0 = 300. Right: Associated 
timescales for shear acceleration, calculated for <fr c = 0.1 [rad]. 
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Fig. 2. — Evolution of the critical comoving proton Lorentz factor 7' as a function of the 
radial coordinate, shown for jj, = 300, £ 77 = 1, (3 S — 0.9, and different jet magnetic field 
configurations. Above 7g q the shear acceleration timescale becomes smaller than the shock 
acceleration timescale and proton acceleration by shear becomes more efficient than shock 
acceleration. 



